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The molecule-based magnet [Ru2(C>2CMe)4]3 [Cr(CN)e] contains two weakly-coupled, interpene- 
trating sublattices in a body-centered cubic structure. Although the field-dependent magnetization 
indicates a metamagnetic transition from an antiferromagnet to a paramagnet, the hysteresis loop 
also exhibits a substantial magnetic remanance and coercive field uncharacteristic of a typical meta- 
magnet. We demonstrate that this material behaves like two giant moments with a weak antiferro- 
magnetic coupling and a large energy barrier between the orientations of each moment. Because the 
sublattice moments only weakly depend on field in the transition region, the magnetic correlation 
length can be directly estimated from the magnetization. 



Due to the importance of weak interactions that are 
commonly neglected in solid-state materials, molecule- 
based magnets exhibit a variety of novel behavior [l|, Q • 
One of the most fascinating and puzzling molecule- 
based magnets containing two interpenetrating lattices 
of [Ru 2 (0 2 CMe) 4 ] 3 [Cr(CN) 6 ] (Me = methyl, CH 3 ) @,i 
(Cr(Ru 2 ) 3 ) exhibits an unusual "wasp-waisted" hystere- 
sis loop that is attributed to the weak antiferromagnetic 
(AF) coupling between its two sublattices [9, [a!. The 
shape of the initial M(H) data plotted in Fig.l p}, with 
a shallow slope followed by a rapid rise, is commonly as- 
sociated with a metamagnetic transition between AF and 
paramagnetic (PM) states Q, which has been observed 
both in solid-state systems like FeCl2 and DyAl garnet 
1C| as well as in several other molecule-based magnets 
11]. However, typical metamagnets do not exhibit the 
nonzero remanent magnetization and coercive field found 
in the hysteresis loop of Cr(Ru 2 ) 3 compounds with inter- 
penetrating lattices. 

We argue that a Cr(Ru 2 ) 3 compound with inter- 
penetrating lattices behaves like two macroscopic AF 
coupled moments with a large energy barrier between 
the orientations of each moment. Quasi-one- or two- 
dimensional materials that magnetically order in three 
dimensions due to a weak coupling between sublattices 
are fairly common, even among molecule-based magnets 
[l2l ] . Much rarer are materials where each of the weakly- 
interacting sublattices is fully ordered in three dimen- 



sions. 

A single lattice of Cr(Ru 2 ) 3 



is sketched in the in- 



set to Fig.l, where every pair of Ru ions bridges two 
[Cr(CN) 6 ] 3 ~ (Cr) ions located at the corners of the cu- 
bic unit cell and separated by a; = 13.4 A. While each 
Cr(III) ion has a spin S = 3/2, each Ru 2 dimer is in 
a II /III mixed- valence state with a total spin 5 = 3/2 
0] . Due to the "paddle-wheel" molecular environment 
produced by the surrounding four Me groups, each Ru 2 
moment experiences a strong easy-plane anisotropy with 
D ss 100 K or 8.6 meV [l| Gj]. The easy plane lies 
perpendicular to the axis joining the Ru 2 sites with the 
neighboring Cr ions. Although a single-lattice Cr(Ru 2 ) 3 
compound has been synthesized [f| , the resulting samples 
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FIG. 1: The initial magnetization (virgin curve) of Cr(Ru 2 ) 3 
with interpenetrating sublattices for temperatures up to 30 
K along with the predicted field dependence. The upper left 
inset shows the predicted ground state of a single sublattice 
for classical spins with large anisotropy. 



are amorphous. A Cr(Ru 2 ) 3 compound with interpene- 
trating lattices contains a second identical lattice inserted 
through the open space of the first lattice resulting in a 
body-centered cubic structure. 

Most of the properties of the interpenetrating-lattice 
compound can be explained by a simple model with 
strong easy-plane anisotropy D on the Ru 2 sites, AF 
intra-sublattice exchange J c between neighboring Cr and 
Ru 2 sites on each sublattice, and a weak AF inter- 
sublattice exchange K c between moments on the two sub- 
lattices. The coupling K c is the sum of the dipolar en- 
ergy and the superexchange interaction through weakly- 
overlapping molecular orbitals on the two sublattices [15| . 

The key to understanding Cr(Ru 2 ) 3 with interpene- 
trating lattices is to construct the correct ground state for 
the single-lattice compound. Each Ru 2 spin pair in Fig.l 
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is labeled as a (along the x axis), b (along y), or c (along 
z). In the classical limit with infinite anisotropy, the a, b, 
or c spins must lie the yz, xz, or xy planes, respectively. 
In the quantum case, the Ru2 spins will have small, but 
nonzero components in the classically-forbidden direc- 
tions. For both classical and quantum spins, AF order 
is frustrated by the easy-plane anisotropy. A similar sit- 
uation arises in cubic pyrochlores like Ho2Ti207, where 
ferromagnetic order is frustrated by local [111] anisotropy 

u. 

Both classical and quantum calculations provide the 
same magnetic ground state. The sum of the Ru2 a, 6, 
and c spins points opposite to the Cr spin along one of 
the four diagonals of the cube. Accounting for the two 
orientations of the moment along each diagonal, there 
are eight domains in zero field. For classical spins with 
infinite anisotropy, the net moment along one of the di- 
agonals is M s i = (VH — 1)5 s» lAbS per Cr ion. In 
both classical and quantum calculations, the Ru2 spins 
approach the cubic diagonal as D/J c is reduced. The 
quantum result M s j w 1.21S 1 for the single-lattice mo- 
ment with D/J c = 5 is lower than the classical result 
because the expectation values of the Ru2 spins are re- 
duced in magnitude from 1.5 to about 1.23. 

Based on quantum mean-field calculations with 
D/J c = 5 and J c = 1.72 meV (see below), the single- 
lattice compound experiences a spin-flop transition at 
H s f w 2.78J c fiB ~ 80 T with the field along one of 
the cubic diagonals. Above H s f, the Cr and total Ru2 
moments cant away from the cubic diagonal and, until 
all the moments become aligned ferromagnetically (FM) 
at an even higher field, the magnetization is no 

longer parallel to H. 

Now consider Cr(Ru 2 ) 3 with interpenetrating lattices. 
Assuming that K c -C J c and [IbH *C J c S, the sublattice 
moments M s ;rij (j — 1 or 2) are rigid with the same 
spin configuration as in the inset to Fig.l but point- 
ing along arbitrary cubic diagonals. Thermal equilib- 
rium between the 64 possible configurations {ni,ri2} is 
established within correlated clusters of size £ containing 
Nci ~ 2(£/a;) 3 Cr spins per cluster. Although £ de- 
pends on field and may peak in the transition region, the 
polycrystalline nature of the Cr(Ru 2 ) 3 samples suggests 
that the average value of £ is most important. With Nqt 
assumed independent of field, the energy is given by 

E = N Cl J2{-^BM sl (nu + n 2i )-U 

i 

+ 3K C S 2 m< • n 2i } , (1) 

where H = if m is the magnetic field and we sum over 
clusters i The intra-sublattice exchange J c only enters 
E through the sublattice moments M s i(T), which vanish 
above T c oc J C S 2 . 

The zero-field AF ground state is shown schemati- 
cally in Fig. 2(b) with ni = — n 2 . Depending on the 
field orientation, the AF state becomes unstable to one 
of three possible spin configurations. Configuration III 




FIG. 2: (a) The T = average moment and non-interacting 
linear susceptibility XnmtH/2 and (b) the fractions of the 
AF state and configurations I, II, and III versus /j,bH/K c S 
assuming classical spins and large anisotropy. 



with ni ■ ri2 = —1/3 appears at the critical field H c = 
V3K c S 2 /(n B M s i) when H is parallel to [111]. In con- 
figuration II, ni • 112 = 1/3 and in the PM state labeled 
configuration I, ni — ri2. The stable phases are shown 
in Fig. 2 along with the fraction of configurations I, II, 
and III when averaged over all field directions m. No- 
tice that the AF configuration vanishes only with the 
appearance of the PM state. In zero field, M av obeys the 
simple functional form M av /M sat = 1 — (H c /H) 2 plotted 
in Fig.2(a). 

We must modify Eq.|T|) to account for the small dis- 
tortion of the single-lattice ground state with field. That 
distortion is responsible for the small linear susceptibil- 
ity observed within the AF state at low temperature and 
small fields, and for the even smaller differential suscepti- 
bility 2nBdM j 'dH observed within the PM state at high 
fields. For classical spins and large anisotropy, the mag- 
netic ground state with = m cannot be deformed 
until very large fields, when the Cr spins tilt away from 
the — m direction. But for rij = — m, the Ru2 spins can 
easily cant towards the field direction. 

So the simplest form for the susceptibility of sublattice 
j is Xj — Xsl s i n2 (#i/2) where 9j = cos~ 1 (n :) - • m) is the 
angle between the field and moment directions. Hence, 
Xj = when 9j = and Xj = Xsl when 8j = tt. The 
non-interacting susceptibility of the magnetic configura- 



3 



7000r 

6000^ 
5000^ 
4000 

3000; 

2000: 
1000 

0: 

o ■-: 
E 5 

^ 4 
E 3 

& z 

0, 



(a) 



(b) 



18000 

16000 

14000 

12000 

10000 

8000 

6000 

2.5x10 
2.0x10 
1.5x10 
1.0x10 
.5x10 



(8 



3 
c 

O 

Ct) 

s 
3 
o 



7s 

3 

3 

3 CD 
< 

3 
3 



5 1 1 5 20 25 3 

T(K) 



FIG. 3: The fitting parameters versus temperature: (a) 
the number iVcr(T) of Cr atoms in each magnetic cluser 
(the dashed curve gives the fit N Cl = N [(T C - T)/T c ]~ Av 
with v = 0.5) and the saturation magnetization 2/j,BM aa t(T) 
and (b) the single- lattice susceptibility Xsi(T) and the inter- 
sublattice interaction K C (T). 



tion {m, n 2 } is then given by Xnint = Xi +X2 per pair of 
Cr atoms and the additional linear term in the magneti- 
zation is NcrXnintH/2. The extra term -N Cr XmntH 2 /4 
must then be added to the energy E of Eq.([T]) for each 
cluster. 

In the AF state with ni = — 112, 62 = #1 + 7r and 
Xmnt = Xsi- In the PM state with 9 1 = 6 2 , Xmm = 
2xsi sin 2 (9\ /2) vanishes when 6% = but is nonzero when 
H points away from a cubic diagonal. Averaged over 
all field directions, the susceptibility in the PM phase 
is Xnint ~ 0.134xsi- So in agreement with experiments 
0) 0] ) the differential susceptibility at high fields is much 
smaller than the the linear susceptibility at low fields. 

Using this model, we evaluate the partition function 
by summing over the 64 configurations of {ni,ri2} and 
the magnetization by averaging over field directions. For 
every temperature, the average magnetization depends 
on four parameters: the single-lattice susceptibility Xsi, 
the sublattice moment M s i , the weak AF interaction K c 
between sublattices, and the number iVo of Cr atoms 
within each cluster, half belonging to each sublattice. 
The fits are provided in Fig.l and the resulting parame- 
ters are plotted in Fig. 3. 

In Fig. 3 (a), Nc r {T) is consistent with the functional 
form N Cr (T) = N a [(T c - T)/T c ]~ 3v , where N » 725, 
T c = 39 K, and v = 0.5, which is the mean-field value 
for the critical exponent. Since the functional form 



£(T) cx [(T c - T)/T c ]-" is only expected close to T c , it 
is remarkable that deviations from this form are only 
apparent at 10 K. Other evidence for the mean-field na- 
ture of the phase transition in Cr(Ru 2 ) 3 comes from the 
magnetization of the single-lattice compound :4| , which 
although amorphous closely follows a Brillouin function 
and vanishes near T c like (T c — T) 1 / 2 . 

The "saturation" magnetization 2fiBM sat (T) = 
\/3/j,bMsi(T) per Cr atom plotted in Fig. 3(a) gives the 
average magnetization when the two nondistorted sub- 
lattice moments are aligned. Nevertheless, the magneti- 
zation continues to rise with increasing field due to the 
deformation of the ferrimagnetic ground state of each 
sublattice. Indeed, Vos et al. [H found that the mag- 
netization at 5 K rises from 16,937 emu Oe/mol at 0.5 
T to 20,800 emu Oe/mol at 5 T, far above the "satu- 
ration" value of 17,286 emu Oe/mol. In rough agree- 
ment with the assumed form for the non-interacting sus- 
ceptibility Xmnt, the high-field differential susceptibility 
of 0.086 emu/mol is about nine times smaller than the 
low-field linear susceptibility of 0.81 emu/mol. Notice 
that the 30 K clusters with Ncr = 6464 have correlation 
lengths £ <~ 15a/ = 200 A and sublattice magnetizations 
(i B N C rM s i(T = 30K)n 4500/i S . 

At T — 5 K, the saturation magnetization is signif- 
icantly lower than the classical result with D = 00. 
But the single-lattice moment M s i for quantum spins de- 
creases rapidly with increasing D/J c from its maximal 
value of 2S per Cr ion when D = 0. The fitted sat- 
uration value can be used to estimate that D/J c w 5 
and J c « 1.72 meV. Another estimate for J c comes from 
the mean-field result T c k, 4.05J c = 33 K for the tran- 
sition temperature with D/J c — 5, yielding the smaller 
value J c 0.70 meV. However, mean-field theory typi- 
cally overestimates T c and underestimates J c . So at low 
temperatures, the inter-sublattice coupling K c is roughly 
600 times smaller than the intra-sublattice coupling J c . 

Plotted in Fig. 3(b), the sublattice susceptibility Xsi 
rises sharply with increasing temperature and is roughly 
proportional to Nc r - The inter-sublattice coupling K C (T) 
falls off almost linearly with temperature below 25 K. 
If K c were short ranged and coupled sublattice mo- 
ments at distinct points in space, we would expect that 
K C (T) ex. M s i(T) 2 . The temperature dependence of 
K c (T)/M s i(T) 2 may be ascribed to the complex and 
long-ranged interaction between the two sublattices, each 
of which contains two species of magnetic ions with dif- 
ferent temperature-dependent average moments. 

Three of the four fitting parameters in this model, 
M s i(T), Ncr(T), and Xsi(T), are properties of the single- 
lattice compound; only K C (T) reflects the presence of two 
sublattices. Above about 30 K, the fits break down pri- 
marily because the linear terms Xsi{T) sin 2 (0/2)7J in the 
sublattice magnetization are no longer small compared 
to the zero-field value of 2/ibM s ;(T). Consequently, the 
inter-sublattice coupling K c effectively depends on field 
as well as on temperature. 

Despite the success of the fits below 30 K, the lower- 
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right inset to Fig.l suggests one limitation of our model. 
Although the correlation length £ must reach a maximum 
near H C (T) = V3K C (T)S 2 /{^ B M s i{T)) as fluctuations 
soften, our model assumes that £(T) is independent of 
field. Consequently, the curves in Fig. 1 lie slightly above 
the data points near H C (T), especially at lower temper- 
atures. A correlation length f (T, H) that peaks in the 
vicinity of H C (T) would produce even better agreement 
with the experimental data. Because the thermal aver- 
ages are assumed to be independent from one cluster to 
another, the model magnetization is a continuous func- 
tion of H for any nonzero temperature and for any orien- 
tation of the magnetic field. Due to the polycrystallinc 
sample, we cannot say whether the experimental mag- 
netization of Cr(Ru 2 ) 3 experiences a jump near H C (T) 
at nonzero temperatures, as observed in single crystals 
of conventional metamagnets 0, 03 up to about 85% of 
the Neel temperature. 

The interpenetrating-lattice compound Cr(Ru 2 ) 3 may 
be the only known material with two or more weakly- 
interacting three-dimensional sublattices. Another 
molecule-based magnet that bears some similarities to 
Cr(Ru 2 ) 3 is methylamine chrome alums (l7| . which con- 
tains Cr(III) ions in a face-centered cubic structure. Al- 
though weak AF dipolar interactions couple the four sub- 
lattice moments, the intra-sublattice exchange is negligi- 
ble in chrome alums. So the three-dimensional ordering 
at about 0.02 K is also produced by the dipolar interac- 
tions between magnetic sublattices. 

Most of the puzzling properties of Cr(Ru 2 ) 3 with in- 
terpenetrating lattices can be explained by two macro- 



scopic moments that are weakly AF coupled with a 
large energy barrier between the orientations of each 
moment. If the individual moments within each sub- 
lattice rotate together, then the energy barrier between 
moments in the [111] and [111] directions for classical 
spins at T = is A(0) = 0. 07 J c S 2 N Cl . At nonzero tem- 
peratures, we expect that A(T) oc N Cl -(T)M s i(T) 2 J c ~ 
N J C [(T C - T)/T c ]- 1 /2. So A(T) will rise very close to 
Tl, which may explain the anomalous AC susceptibility 
[3| observed in these materials. 

Our work provides several new predictions for 
Cr(Ru 2 ) 3 with interpenetrating lattices. In the magnetic 
ground state, the total moment of each sublattice must 
lie along one of the four cubic diagonals. This can be ver- 
ified by neutron scattering, either on a deuterated sample 
or with a large enough polycrystal. With a single crystal, 
the transitions between the AF and configurations I, II, 
or III should be observable for different orientations of 
the field. The predicted correlation length £(T) can be 
verified independently by fitting the elastic peaks mea- 
sured with neutron scattering, even on a polycrystalline 
sample. We hope that this paper will inspire future work 
on this fascinating system, including measurements on 
single crystals when they become available. 
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